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Abstract 

The latest top quark mass measurements by the CDF and D0 experiments are 
presented here. The mass has been determined in the dilepton (tt — *■ e/i, ee, /x/^ 
+ jets and lepton plus jets (tt — + e or /i + jets final states. The 

most accurate single result from lepton plus jets channel is 173.5l3'g(stat. + 
Jet Energy Scale Systematic)±1.3(syst.) GeV/c^, which is better than the 
combined CDF and D0 Run I average. A preliminary and unofficial average 
of the best experimental Run II results gives Mtop = 172.7 ± 3.5 GeV/c^. 



1 Introduction 



Since the first evidence in 19941^ and the discovery of the top quark in 
, the CDF and D0 CoUaborations invested a lot of work to determine 
its properties, specially the value of its mass, which is a fundamental parameter 
of the Standard Model (SM). The ongoing Run II, with the upgraded Fermilab 
Tevatron collider and CDF and D0 detectors, gives new possibilities for a 
precise measurement of the top mass. Due to its large mass, corresponding to 
a Yukawa coupling of order unity, one may suspect that the top quark may 
have a special role in the elcctroweak symmetry breaking. In addition, due to 
its significant contribution to high order radiative corrections of a number of 
electroweak observables, a precise measurement of the top quark mass provides 
a tighter constrain on the Higgs mass®. 

This paper reports on the latest CDF and D0 top quark mass results 
which are based on about 318 pb^^ (CDF) and 219 pb"^ (D0 ) of data from 
the first two years of the Tevatron Run II (2002 to 2004). Another paper, pre- 
sented on this conference ^ , summarized the top quark kinematics properties 
including its recently measured production cross section at the center of mass 
energy of v^=1.96 TeV, 

At this Tevatron energy, top quarks are produced generally in pairs from 
the processes qq — > tt (in ~85% of the cases) and gg tt (in ~15% of the 
cases). Top can be produced as a single quark by electroweak interactions, 

by W-gluon fusion or virtual W* production in the s-channelS', but with a 
smaller cross section. At this time, no signal has been observed from single top 
processes and they are not expected to be utilized for a precise mass measure- 
ment. 

In the Standard Model the branching ratio of the decay t bW is nearly 
100%. When a tt pair is produced, each of the W-bosons can decay into either 
a charged lepton and a neutrino (branching ratio of 1/9 for each lepton family) 
or into a qq quark pair (branching ratio of 2/3). This allows us to classify the 
final states as: 

• Dilepton final state, when both W's from the tt pair decay leptonically. 
This state is characterized by two high Pt charged leptons, two jets 
from bb quarks ^ and significant missing transverse energy (^t) from the 
neutrinos. 

• Lepton plus jets final state, when one W boson decays leptonically and 
the other one hadronically. This state contains one high Pt charged 

"'^ Errors in jet reconstruction and gluon radiation in the event may make the 
observed number of jets smaller or larger. This statement is valid for all final 
states. 
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Figure 1: The summary of the top mass "evaluations" and direct measurements 
versus time. See the text for explanation of the points. 

lepton, four jets and significant EIt- 

• All hadronic final state, when both W's decay hadronically. This state is 
characterized by 8 jets, two of which arc from 6-quarks. 

In Run I, CDF and D0 used all of these signatures for the top mass 
measurements. At this time, the Run II top mass analyses from the all-hadronic 
channel are still in progress and will not be reported in this paper. 

Figure n shows how our knowledge on the top quark mass improved with 
time. The diamonds represent the indirect determinations from fits to the 
electroweak observables ^3 . The curves in the upper left corner of the figure 
represent the limits from direct searches in e^e~ and pp machines. The Run I 
CDF and D0 results are presented with squares. The filled circles are the 
new Tevatron Run II results. The band represents the average Run I top mass 
results. One observes that the recent Tevatron results have equal or better 
accuracy than the Run I world average. 
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2 Top kinematics and mass reconstruction methods 



The kinematics of the events in the lepton plus jets final state is over- 
constrained. In this channel, the number of measured quantities and the num- 
ber of applicable energy-conservation equations, from five production and decay 
vertices, are larger than the number of non-measured kinematical event param- 
eters. This feature allows for a complete reconstruction of the four-momenta 
of the final state particles in the event, for example using the two constrain 
kinematical fit (2CF), and for a reconstruction of the top mass event- by- event. 
In this type of analysis there is an ambiguity in how to assign the four leading 
jets to the two b quarks and two light quarks coming from the tt system. If 
none of these jets is tagged by b taggers ^, there are 12 different ways of as- 
signing jets to the 4 partons. Combining with the ambiguity from solving of 
a quadratic equation on P!f, there are 24 different values of mtop returned by 
the fit. The combinations are reduced to 12 or to 4 if one or two of the jets are 
selected from b-taggers. 

The variety of the CDF and D0 analyses in the lepton plus jet sam- 
ple is relatively large. However, all analyses can be separated in three major 
categories: 

1. Template type analyses, where the reconstructed mass distribution from 
the data is compared with expected distributions from Monte Carlo gener- 
ated signal (mass-dependent) and background. In these types of analyses, 
all events are weighted equally. By doing so one neglects the additional 
information coming from a different mass resolution in single events. In 
addition no use is made of possible information from the dynamic of the 
process which can be assumed to be known. Typical examples for this 
type of analysis are the CDF and D0 Run iCSini) 

and the latest, most 

accurate CDF Run II lepton plus jets analysis, described in Section 3.1. 

2. Matrix Element type analyses, originally proposed by Dalitz and Gold- 
stein and independently by Kondo^^. These methods calculate the 
posterior probability, given the known production cross section, for every 
event with measured kinematic properties, to originate from a pp — > tt 
process. A typical example is the D0 Matrix Element analysis which 

is the base of the best Run I top mass measurement . In Run II 

CDF uses a method, proposed by K.Kondo , called Dynamic Like- 
lihood Analysis, which differs in the way that the normalization of the 
differential cross section is performed. 

^CDF and D0 use taggers based on either displaced vertices (Secondary 
Vertex Tagging, SVX, for example see or on low Pt electrons or muons 
from the 6-quark semileptonic decays (Soft Lepton Tagging, SLT^^) 



3. A mixture between methods 1) and 2). For example the tt event is recon- 
structed using the kinematic algorithms similar to the template analyses 
but an event-by-event probability of each kinematic reconstruction is ex- 
ploited as a weight (for example exp{—^)). A typical example of this 

category is the D0 Ideogram Analysis 1^ . 

3 Lepton plus jet channel 

3.1 CDF result 

For the time being the most accurate top quark mass measurement comes 
from the lepton plus jets channel. This channel combines the benefits of good 
signal to background ratio, the possibility to reconstruct the top quark mass 
event by event with a relatively small combinatorial effect, and a large branch- 
ing fraction. In brief, we discuss below the main selection criteria for this 
channel. 

Lepton plus jets events have the signature pp —^ tiX —^ ivhqq bX. The 
characteristics of this final state begins with the identification of one isolated 
central high energy lepton (e with Et > 20 GeV or /i with Pt >20 GeV) and 
|,7|<l3. ^ 

Assuming that the lepton is coming from W boson decay, a companion 
neutrino should exist. This would spoil the balance of the observed momentum 
in the transverse plane. The missing transverse energy (E!t) is constructed by 
adding the calorimeter energy vectors in the plane transverse to the beam. The 
calorimeter clusters identified as jets are corrected for detector response and for 
multiple pp interactions. In muon events the is computed using the muon 
momentum measured by the track instead of the muon calorimeter signal. 

In order to fully reconstruct the tt system, at least four central jets |?7| < 2 
are required in the system. The SVX tagging algorithm is run over the lead- 
ing jets {Et >15 GeV): some of them may be identified as b-jets. To obtain 
maximum statistical benefit from the event sample it is helpful to decompose 
it into several classes of events which are expected to have different signal-to- 
background ratios and top mass resolutions. The Monte Carlo studies showed 
that an optimal partitioning is obtained splitting the sample into four statisti- 
cally independent categories: events with double SVX tags (2SVX), events with 
single SVX tag and tight forth jet cut (four jets with Et >15 GeV, ISXVT), 
events with single SVX tag and loose forth jet cut (15 > Et >8 GeV for the 4*'' 
jet, ISXVL) and finally events without tags (0-tag). Since the last sample has 
a high background contamination compared to SVX tagged ones, an additional 



^A complete description of the lepton selection, including all cuts used, can 
be found elsewhere^!' 



Table 1: In the first four columns from left to right: lepton plus jets suhsamples 
used in the top quark mass analysis, number of events in each sample, S/B 
ratio, and a summary of the jet energy cut selection are presented. A total of 
165 tt candidates were selected. The last column summarizes the background 
fractions in % from W + light quark, Whh + Wcc + Wc and QCD multijet 
events (left to right). 



Data 
Subsample 


Number 
of events 


S/B 


Jet Et cuts (GeV) 
jets 1-3 (4*/i jet ) 


Bckg. type and 
fraction in % 


2SVX 


25 


10.6/1 


Et >15 ( Et >15) 


21/59/10 


ISVXT 


63 


3.7/1 


Et >15 ( Et >15) 


17/38/22 


ISVXL 


33 


1.1/1 


Et >15 ( 15 > >8) 


29/48/14 


0-tag 


44 


0.9/1 


Et >21 ( Et >21) 


75/3/20 



optimization of the jet Et cuts ( Et > 21 GeV) was performed. A total of 165 
tt candidates were selected from 318 pb^^ of data. 

The dominant backgrounds in all samples are direct W plus multijet pro- 
duction, including heavy flavour production, and QCD multijet events where 
one jet is misidentified as a lepton. Additional small backgrounds are due 
to WW/WZ and single top production. The amount and composition of the 
background depends on the sample. In the case of the 2SVX sample, the 
W^bb -I- VFcc -I- VFc background dominates (~60%) while in the case of the 0-tag 
sample, the W plus light quark production is responsible for ^75% of back- 
ground. The information for these four subsamples, including the dominant 
type and background fraction, is summarized in Tabled 

Each event, either from data or MC samples, is fitted to the hypothesis 
tiX liybqq bX . We use four kinematic constraints, as a consequence of the 
assumed lepton plus jets event structure {Mi^, = Mjj = M\y and Mii,h = Mjjb 
= Mtop). The fitting procedure runs over all possible 24 combinations of as- 
signing the four leading jets to the 6, b and W ^ qq partons (the order of the 
pair assigned to the W is irrelevant). If one or two of the four leading jets are 
tagged as a 6-jets, they are assigned to the 6-partons and the number of ex- 
plored combinations is correspondingly smaller. All solutions with < 9 (cut 
optimized on the MC studies) are accepted. The top mass value corresponding 
to the combination with the minimum is picked as the mass value indicated 
by the event. 

The events from the MC samples are used to produce probability density 
distributions or so-called templates. In case of the signal MC, these distribu- 
tions are parameterized as a function of reconstructed and input top masses. 
On the other hand, the background probability density distributions are pa- 
rameterized only as a function of the reconstructed mass. The likelihood of 




Figure 2: Left: reconstructed top mass p.d.fs for input top masses from 14-5 to 
205 GeV/c^ and JES — 0. Right: reconstructed mass p.d.fs of the dijet system 
attributed to the W for different JES values in the range of —3(7 — + 3(t and 
input top mass of 180 GeV/c^ . 



each subsample uses the parameterized signal and background probabiUty den- 
sity functions (p.d.f) to evaluate the dependence of the likelihood on the input 
top mass. 

To reduce the dominant systematic error coming from jet energy scale 
(JES) the latest CDF template analysis exploits the fact that the global 
JES scale can be determined from the decay W ^ qq . MC studies have shown 
that this technique provides a 22% reduction in the JES uncertainty. Similar 
template distributions as for the kinematically reconstructed top mass are built 
for the dijet mass, with the exceptions of removing the < 9 requirement, 
exploiting all the possible jet-to-parton assignments in the event. Examples of 
top and dijet reconstructed masses p.d.fs are shown on Figure 2 left and right 
respectively. 

The reconstructed top and W dijet mass values for every data event are 
simultaneously compared to the p.d.fs from signal and background sources 



Table 2: The systematic uncertainties in the CDF lepton plus jets top quark 
mass measurement. 



Source 


(GeV/c2) 


b-jets modeling 


0.6 


Method 


0.5 


Initial state radiation 


0.4 


Final state radiation 


0.6 


Shape of background spectrum 


0.8 


6-tag bias 


0.1 


Parton distribution functions 


0.3 


Monte Carlo generators 


0.2 


MC statistics 


0.3 


Total 


1.5 



performing an unbinned likelihood fit. The fit finds a maximum likelihood 
value according to: the expected numbers of signal and background events in 
each subsample, the JES and the true top quark mass {Mtop)- Only Mtop is 
a free parameter in the likelihood fit, the other are constrained within their 
uncertainties. For each subsample, the likelihood has the following form: 

^ ~ ^ shape ■ ^ shape ' ^counting ' i-bg- j 

In JQ, the main information on the top quark mass is hidden in the 
term C^f^"^^. It gives the probability for an event with reconstructed top mass 
77i[qp to come from true top mass Mtop- AH other terms constrain the JES 
(^Thape)' ^^'^ number of observed events {Ccounting) and number of expected 
backgrounds in the subsample (Cbg) and help to reduce the statistical and 
systematical uncertainties returned by the fit. 

The systematic uncertainties are summarized in Table |21 For each sys- 
tematic source, the relevant parameters are varied by ±la in the tt MC sample 
with Mtop=l78 GeV/c^ and sets of fake events are generated. These fake events 
are reconstructed in the same way as normal events. This procedure is called 
"pseudoexperiments" (PE). It propagates the ilcr effects to a shift in the top 
mass relative to the result from the nominal sample. 

The reconstructed top masses in the four subsamples with overlaid best fit 
for the signal and background MC expectation are shown in Figure 3. The com- 
bined fit for all lepton plus jets events returned Mtop — 173.5^3 g(stat.-|-JES) 
±1.3(syst.) GeV/c^ and JES = -0.10to;gi(stat.+syst.). This is the most precise 
single measurement available to date, better than the average Run I result. 
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Figure 3: The reconstructed mass distribution (histogram) for each lepton 
plus jets CDF suhsamples is overlaid with the result of the likelihood fit (sig- 
nal+background, hatched area). The cross hatched area represents only the 
background. 



3.2 D0 Result 

D0 has measured the top quark mass in the lepton plus jets channel as 
well. The utilized data set corresponds to an integrated luminosity of approx- 
imately 229 pb"^ collected between April 2002 and March 2004. 

The event selection criteria are similar to those used in CDF. As a first 
selection step, an identification of the a high Pt isolated electron or muon ac- 
companied by substantial large > 20 GeV is required. The isolated electron 
(muon) candidate should have Pt >20 GeV, satisfy a pseudo-rapidity cut of 
I77I <1.1 {\ri\ <2.0) and tight quality conditions. These initial selections provide 
the data sample. 

Two separated analyses, b-tagged and topological, are performed on this 
sample. In the b-tagged analysis, to reconstruct the top mass the events are 
additionally selected to have at least 4 jets with Pt > 15 GeV and |ry| <2.5. A 



further requirement of identification of one or more jets as 6-jets is made. A 
jet is b-tagged based on the reconstruction of the secondary verticies using the 
charged particle tracks associated with it. 49 (29) e+jets (/i + jets) b-tagged 
events survive all cuts and are kinematically fitted to the tt hypothesis. In 
42 (27) electron (muon) events the kinematic fit converged in a configuration 
where the lowest solution is consistent with b-tagged jet permutation. 

In the second analysis the information of the b-tagger is not exploited. 
To increase the signal to background ratio several modifications of the selection 
cuts are applied. For example the transverse momenta of the first four jets are 
increased to 20 GeV. There are 87 e-t-jets and 80 /i-|-jets events left after this 
requirement. Next, using the specific kinematics of the tt events, a discriminant 
(D) was constructed. It is designed to use variables which are uncorrelated or 

minimally correlated with the top quark mass Four topological variables 
are considered: 

• - missing transverse energy which comes from the neutrino of the W 
leptonic decay. 

• A - aplanarity of the event. It exploits the fact that the decay products 
from a massive particle have large aplanarity. 

• Hrp2 - the ratio of the scalar sum of the Pt of the jets, excluding the 
leading jet, and the scalar sum of \pz\ of the jets, the lepton and of the 
reconstructed neutrino. 

• Krp^^^ - a measure of the jet separation folded with the Et of the recon- 
structed leptonic W boson. 

Figure 4 shows the simulated discriminant distribution for signal and 
background. A cut of D > 0.4 is imposed to select the signal rich region. After 
the kinematical fit at least one jet permutation is required to have < 10. 

Similar to the CDF lepton plus jets analysis, two dominant sources of 
background are accounted for: W plus multijet production, including heavy 
flavour, and QCD multijet events where one of the jets is misidentified as 
a lepton and there is significant E^t imbalance in the event due to detector 
resolution. 

The systematic uncertainties of both analyses are summarized in Table|3| 
The main contributions are due to JES, gluon radiation (initial state-ISR and 
final state-FSR), and the MC tt signal modeling. 

As expected, the dominant systematic uncertainty originates from the 
JES. The impact of JES on the reconstructed top mass was evaluated by scaling 
the jet energies by ±la for both signal and background in the MC simulation. 
The uncertainty on the JES was conservatively assumed to be 5% for the Ei^^ > 
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Figure 4: Discriminant D for tt (solid dark line) and background (solid light 
line) events from MC simulation. 

30 GeV. For jets with E^* < 30 GeV, the JES uncertainty decreases hnearly 
as cr = 30% - 25% x (£;f^V30) GeV. 

Next in importance to the JES is the systematic uncertainty coming from 
gluon radiation. Regardless of which jet permutation is used, the fitted mass 
will not be correct if the a radiated gluon is one in the four leading jets in the 
event. To understand how this affects the tt reconstruction, MC events with 
only four partons hadronizing and forming four jets were compared to events 
where one of the leading 4 jets comes from gluon radiation. A small deviation of 
~ 0.2 GeV/c^ from the nominal top mass is observed when the events without 
gluon radiation are reconstructed. The difference becomes ^ 2.4 GeV/c^ when 
one of the leading jets is a radiated gluon. 

In this analyses the model of the kinematic properties of the events is 
taken directly from MC simulation. Therefore some deficiencies in the MC 
model may lead to a substantial bias in the mass reconstruction. In order 
to perform a conservative estimate of this possible effect, in addition to the 
nominal sample for tt signal a complementary sample was generated where an 




Table 3: Systematic uncertainties in the D0 lepton plus jets top quark mass 
measurement. The uncertainty on JES, Gluon Radiation and Signal Model are 
the dominant sources of error on the mass. 



Source 




AMtop 








(GeV/c2) 






b-taggcd analy 


3is 


Topological Analysis 


JES 


+4.7/-5.3 






+6.86.5 


Jet Resolution 


±0.9 






±0.9 


Gluon Radiation 


±2.4 






±2.6 


Signal Model 


+2.3 






+2.3 


Background Model 


± 0.8 






± 0.7 


6-tagging 


±0.7 






N/A 


Calibration 


±0.5 






± 0.5 


Trigger bias 


±0.5 






± 0.5 


MC statistics 


±0.5 






± 0.5 


Total 


±6.0 


+7.8-7.1 



additional parton is produced in association with the tt pair. The cross section 
for this process is approximately two times smaller than the cross section for 
the tt production. By analyzing this sample, an uncertainty of +2.3 GcV/c^ 
due to the uncertainty on signal modeling is assigned to the analyses. All other 
possible systematic effects turned out to be relatively small, at the level of 
0.5-0.7 GeV/c2. 

The distributions of the fitted masses and -ln(£) curves are shown in 
Figure 5. The top two figures show the result from the ^-tagged analysis while 
the bottom two represent the topological one. Taking into account the output 
from the binned likelihood fit and the systematic uncertainties the final result 
for the analyses is Mtop = 170.6 ±4.2(stat.) ±6.0(syst.) GeV/c^ (^-tagged 
analysis) and Mtop = 169.9±5.8(stat.)l7;^(syst.) GeV/c^ for the topological 
one. 

4 Dilepton Channel 

4.1 CDF result 

CDF has several independent dilepton analyses which are found to return 
consistent values for the top mass. Since this sample has good signal to back- 
ground ratio (~ 4/1) one is stimulated to invent ingenious ways to reconstruct 
the events and extract Mtop- 

The event selection criteria are similar as in the lepton plus jets channel. 
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Figure 5: The result from the binned likelihood fit for the b-tagged (topological) 
analysis is presented on the upper (lower) left plot. The dots represent the data, 
the solid line is the fitted ti plus background and the dashed line is background 
only, normalized to the fraction returned by the fit. The right plots show the 
-ln(C) curves. 



Two or more central jets with Et > 15 GeV are required. A loose criterion is 
applied to the second lepton - it must have opposite charge but isolation is not 
mandatory. For the missing transverse energy the cut is increased to >25 
GeV since two neutrinos are supposed to be presented in the event. If EJt <50 
GeV, a requirement for the angle between and the nearest lepton or jet to 
be A(j) >20° is imposed. Also the transverse energy sum, Ht, has to be more 
than 200 GeV. Events due to cosmic rays, conversions or Z bosons are rejected. 

Four major backgrounds are taken into account: di-boson plus jet pro- 
duction, W plus jets where one of the jet is faking a lepton, and Drell-Yan 
production, specially Z/7 — > rr. 33 events passed all cuts with an expected 
background of 11.6±2.1 events. 

In contrast to the lepton plus jets mode, in the dilepton case due to 
the presence of two neutrinos the kinematics is not constrained. The number 



of non-measured kinematical variables is larger by one than the number of 
kinematic constraints {—ICF). Obviously, it is impossible to single out only 
one solution per event. We may take some event parameter (i?) as known 
in order to constrain the kinematics and then vary R to determine a set of 
solutions. In order to determine a preferred mass, every solution should have 
a weight attached to it. 

The minimal requirement in the case of —ICF kinematics is to use a two 
dimensional vector as R. We chose the azimuthal angles of the two neutrino 
momenta R ~ {(t>ui, 0i/2) and create a net of solutions in the (^i/i, </>iy2) plane. 

For every point of the {(t)ui,4'u2) plane we have 8 solutions. Two of them 
are generated by the two possibilities of associating the two charged leptons 
to the two leading jets which are assumed to originate from the bb partons. 
The four other solutions are generated from the four ways of associating each 
neutrino to two Pz momenta, satisfying the W decay kinematics. We select the 
minimal solution for every point of the net for further use in our analysis. 

Using the value from a minimization we weight the selected solutions 

2 

by e(— ^). This is done in order to suppress the solutions which have worse 
compliance with the fit hypothesis. 

The final extraction of the top quark mass from a sample of dilepton 
candidates is provided by the unbinned likelihood fit. The expected signal 
and background p.d.fs are obtained using Monte Carlo samples with detector 
simulation. The background-constrained fit (N;,=11.6±2.1) returns: Mtop = 
169.8 ±11 GeV/c^ with 23.4±y signal events. The left plot in Fig.Elshows 
the fitted mass distribution. The insert shows the mass dependence of the 
negative log-likelihood function. The right plot represents the error distribution 
for Monte-Carlo simulated experiments, where the arrows indicate the data 
result. 

We also performed a fit without constraining the number of background 
events. This fit returns Mtop = 169.2±|| GeV/c^, with 33.0±i:g signal events 
and 0.0±Q Q background events. 

4.2 D0 result 

To reconstruct the top quark mass in the dilepton channel, D0 follows 
the ideas proposed by Dalitz and Goldstein in ^'^K The analysis uses about 
230 pb~^ of data. The initial selection includes: 

• two leptons, electron or muon, with p^ > 15 GeV in the pseudorapidity 
regions |?7(e)| <1.1 or 1.5< \r]{e)\ <2.5 for the electron and \ri{fi)\ <2.0 
for the muon. In e-/i events a separation cut of •^/A^^TT'Att^ >0.25 is 
applied; 

• two or more jets with Pt > 20 GeV in the pseudorapidity region \r]^'^^ \ <2.5; 




Reconstructed mass (GeV/c^) Statistical error (GeV/c^) 



Figure 6: Left: two-component constrained fit to the dilepton sample. The 
cross hatched area corresponds to the background returned by the fit and the 
line-shaded area is the sum of background and signal. The insert shows the 
mass- dependent negative log-likelihood used in the fit. Right: positive and neg- 
ative error distributions returned by the fit in pseudo experiments. The arrows 
indicate the errors from the data fit. 

• large missing >25 GeV. However the EIt cut is varied for di-electron 
or di-muon events depending on the ee or /i/i invariant mass; 

• veto on the Z — > ee, ^/z events; 

• a cut A0(^, ESt) > 0.25 rejects the events where the ESt and ^ vectors are 
close to each other in the transverse plane; 

• Ht > 140 GeV, where Ht is the scalar transverse momentum sum of the 
larger of the two lepton pxS and of all jets over 15 GeV. 

• for ee events a additional sphericity > 0.15 cut is applied. 

8 e^, 5 ee and fifi events satisfy all requirements, when 6.2±0.6, 2.8±0.3 
and 2.9±0.6 events are correspondingly expected. 



The D0 analysis method can be summarized as foUows. The momenta of 
the two highest Pt jets in the event are assigned to the bb from the decay of tt 
quarks. Then a Ukehhood to hypothesized values of the top mass in the region 
of 80^280 GeV is determined. For each event a solution is found when the 
pairs of tt momenta are consistent with the observed lepton and jet momenta 
and ESt- a weight to each solution is assigned as 

W = f{x)f{x)piEnMtop)p{EI\Mtop), (2) 

where f{x) (/(a;)) is the parton distribution function for the proton (anti- 
proton) and the initial quark (anti-quark) is carrying a momentum fraction x 
(x). p(E^\Mtop) denotes the probability for the top (anti-top) quark with a 
mass Mt to generate a lepton £ {£) with the observed energy in the top quark 
rest frame. 

There are two ways to assign the two jets to the b and b quarks. In addi- 
tion, for each jet-to-parton assignment, there might be up to four solutions for 
each hypothesized value of the mass, coming from the fact that every neutrino 
may have up to two real solutions for its Pz, satisfying the kinematics. Then 
the likelihood for each value of the top quark mass Mtop is given by the sum 
of the weights Wij over all possible solutions: 

W(MO= E E (3) 

p^ solutions jet assignment 

Up to now there was an implicit assumption that all momenta are mea- 
sured perfectly. Therefore the weight in Q is zero if no exact solutions are 
found. To account for detector resolution the weight calculations are repeated 
with input values for the particle momenta drawn from Gaussian distributions 
with means equal of the measured values and widths corresponding to the detec- 
tor resolution. In addition the '^t, value is recalculated from generated particle 
momenta and a random noise from a normal distribution with mean GeV and 
rms 8 GeV is added. Figure [T] up (down) shows the weight curves for e/i (ee) 
events before (solid line) and after (shadow area) resolution smearing. For each 
event the value of the top quark mass at which the weight curve reaches its 
maximum is used as the estimator of the mass. After that, to extract the most 
probable top mass value from the data sample, a standard template method 
which exploits a binned maximum likelihood fit is applied. The likelihood fit 
returns Mtop = 155. t\t(.stat.) ± 7.{syst.) GeV/c^. The JES uncertainty (5.6 
GeV/c^) dominates the systematic error. 
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Figure 7: Weight curves from 8 e/i events (top) and 5 ee events (bottom). The 
shadow histograms show the weight curves with resolution smearing while the 
open histograms represent the weight curves without resolution smearing. 



5 Summary of the Top Quark Mass Measurements and Run II 
Prospects 



Combining the presently available most accurate Run II CDF and D0 mea- 
surements in the dilepton and Icpton plus jets decay topologies, one finds 
172. 7± 3.5 GeV/c^. This result is unofficial. The average is made by the author 
assuming simple correlations (0 or 1) between the systematic uncertainties in 
the CDF and D0 measurements. 

The expected CDF uncertainty for JES systematics as a function of inte- 
grated luminosity is shown in Figure |S1 left. The right plots shows the total top 
mass error versus integrated luminosity for the CDF lepton plus jet analysis. 
One may conclude that with a Run II integrated luminosity of 8 fb~^ the top 
quark mass could be measured by CDF with a precision of ~2.0 GeV/c^. This 
optimistic forecast is based on the present understanding that both the statisti- 
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Figure 8: The expected JES uncertainty from W—^jj as a function of integrated 
luminosity is shown on the left plot. On the right, the total expected top mass 
uncertainty, from CDF lepton plus jets events as a function of integrated lumi- 
nosity, is shown. 



cal and JES systematic uncertainties will decrease as expected with increasing 
integrated luminosity. 

6 Conclusion 

The top quark CDF resuhs from the Tevatron 2002-2004 Run II, with an 
integrated luminosity of 318 and 230 pb^^ for CDF and D0 are presented. 
The best, up to date, measurement of the top quark mass from the CDF lepton 
plus jets analysis is 173.5^3 g GeV/c^. Combining the CDF and D0 dilcpton 
and lepton plus jets Run II results, the author's average of the top quark mass 
is 172.7 ±3.5 GeV/c^. 
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